Integer Quantum Hall Effect in a Lattice Model Revisited: Kubo Formalism 
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We investigate numerically the integer quantum Hall effect (IQHE) in a two-dimensional square 
lattice with non-interacting electrons in presence of disorder and subjected to uniform magnetic 
field in a direction perpendicular to the lattice plane. We employ nearest-neighbor tight-binding 
Hamiltonian to describe the system, and obtain the longitudinal and transverse conductivities using 
Kubo formalism. The interplay between the magnetic field and disorder is also discussed. Our 
analysis may be helpful in studying IQHE in any discrete lattice model. 
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I. INTRODUCTION 

The discovery of quantum Hall effect in two- 
dimensional (2D) electron systems^ exposed to strong 
perpendicular magnetic field was a triumph of experi- 
mental physics. Immediately after this discovery scien- 
tists from various disciplines were launched into a frenzy 
of activity to understand the underlying physics and also 
to explore its technological importance in designing dif- 
ferent electronic devices. Such efforts have led to the en- 
dowment of a new metrological standard, the resistance 
quantum, h/e 2 , containing two fundamental constants, 
the electronic charge e and the Planck's constant /A The 
scaling theory of localization^ suggests that in absence 
of magnetic field all the states of a 2D disordered non- 
interacting electron system are localized due to quantum 
interference. The time reversal symmetry is broken in 
the presence of magnetic field and a series of Landau 
bands appear due to disorder. There are numerous stud- 
ies in the literature on the localisation problem of the 
Landau levels^, specifically studies based on percolation 
theory! and the calculation of Thouless number—. Within 
each Landau band a central region of extended states is 
flanked on both sides by regions of localized states. For 
a disordered 2D electron gas it is now well-understood 
that central region of extended states play a significant 
role towards the quantization of Hall conductance. In 
the limit of strong disorder or weak magnetic field ex- 
tended states in the Landau bands float up in energyiLIii, 
and, a systematic float-up results a transition from inte- 
ger quantum Hall effect states to insulator oneii. 

The interplay between magnetic field and disorder is 
an essential issue for the phenomenon of IQHE. Breaking 
of translational invariance by impurity potential carries 
an essential need for the quantization of electronic con- 
ductivity in presence of magnetic field since in absence 
of disorder no plateau-like structure appears in the Hall 
conductivity, and the classical Hall result is preserved. 
Appearance of Hall plateaux is independent of physical 
dimension of the sample and this universality of the quan- 
tum Hall phenomenon helps the technological progress of 
semiconductor physics using 2D electron gas. Not only 



that, the shape of the system is irrelevant in this con- 
text. With tremendous experimental successes, the phe- 
nomenon of IQHE has drawn a lot of interest among the- 
oreticians too. Most of the theoretical studies available 



FIG. 1: (Color online). Schematic diagram of a square lattice 
subjected to a perpendicular magnetic field B. 

in the literature are based on the continuum model, ig- 
noring the underlying lattice structure of the sample, to 
explain this phenomenon of IQHE. Although there are a 
few attempts based on discrete lattice modeler—, but a 
comprehensive study of this phenomenon with the frame- 
work of the lattice model is still lacking. In the present 
article we have made an attempt in this direction and 
study in details the behavior of IQHE in a square lattice 
in presence of disorder and address several important is- 
sues. 

With this introductory remarks we organize this paper 
as follows. In section II we describe the model and the- 
oretical technique for the calculation of longitudinal and 
transverse conductivities based on Kubo formalism. The 
numerical results and discussion are included in section 
III. Finally, we draw our conclusions in section IV. 



II. MODEL AND THEORETICAL 
FORMULATION 

A. The Model 

We start with a 2D square lattice (Fig. [T]) which con- 
tains N x and N y number of atomic sites along the x- and 
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y-dircctions, respectively and subject it to a perpendic- 
ular magnetic field B. We employ a tight-binding (TB) 
Hamiltonian to describe the system, and, under nearest- 
neighbor hopping approximation it reads as, 
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where, e m .„ denotes the site energy of an electron at site 
(m,n), m and n being the x and ?/ co-ordinates of the site, 
respectively (setting lattice constant a = 1). In order to 
incorporate impurities in the sample, we choose site en- 
ergies e m ^ n randomly from a 'Box' distribution function 
of width W. Here, t is the hopping integral between 
the neighboring sites either along x- or y-direction and 
c m n ( c m,n) is the creation (annihilation) operator of an 
electron at the site (m,n). In the presence of magnetic 
field B, a phase factor 9 m is introduced into the above 



Hamiltonian (Eq. [T]) and for a particular choice of gauge 
A = (0, Bx), the so-called Landau gauge, it becomes non- 
zero only when an electron moves along the y direction. 
Then this phase factor can be expressed as 6 m = 27rm0, 
where <p is the magnetic flux per plaquette measured in 
units of the elementary flux-quantum 4>o = ch/e. We 
set <j> = l/Q, where Q is an integer and the choice of Q 
should be such that it is commensurate with N y . This re- 
duces the boundary conditions to the traditional periodic 
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B. Linear Response Kubo Formalism 

To obtain the longitudinal and transverse conductivi- 
ties, we use the Kubo formalism which is briefly outlined 
below. The general expression for electrical conductivity 
is written in the form, 
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where, r\ — > + . Here, the indices k and I can be x or 
y. For k = I = x, we get a xx , the so-called longitu- 
dinal conductivity while for the other case we have the 
transverse conductivity, a xy . The states \a) and |/3) are 
the eigenstates of the Hamiltonian (Eq. [1} correspond- 
ing to the energy eigenvalues £ a and £p, respectively, 
and N = N x x N y represents the size of the sample. 
iik is the velocity operator along fc-th (x or y) direc- 
tion and f a (p) = 1/[1 + Exp{(£ a (p) - E F )/k B T)}\ is 
the Fermi distribution function at absolute temperature 
T with Fermi energy Ep. 

When k and I are identical to each other, the fac- 
tor {a\uk\P)(P\ui\a) in Eq. [2] becomes a real and pos- 
itive one, while it becomes purely imaginary when k 
and I arc different. Since the conductivity is a real and 
positive quantity, the first term on the right hand side 
of Eq. [2] gives transverse conductivity {<J xy ), while the 
second term represents longitudinal conductivity (a xx ). 
Separating the transverse and longitudinal parts we can 
write, 
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• Matrix elements of the Velocity operators: 

To calculate the matrix elements of the velocity operators 
ii k w e start with the basic relation^, 



Uk = T7~ [u k ,H] 

in 



(5) 



where, Uk is the displacement operator along A:-th (x or 
y) direction and H is the Hamiltonian operator described 
in Eq. [TJ Expanding and simplifying the above relation 
(Eq. [5]) we get the velocity operators in the following 
second quantized form along the x- and y-directions, 
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Therefore, the textcolorredmatrix elements of the veloc- 
ity operators with respect to the eigenvectors \a) and \0) 
become, 
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where, the eigenvectors look like, 
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Here, |p,g)'s are the Wannier states and q and a^ 9 's 
are the corresponding coefficients. Taking the contribu- 
tions from all these states we find the matrix elements of 
the velocity operators along the x- and y-directions. 

Substituting the above expressions (Eqs. [5]and[n]) into 
Eqs. [3] and [4] we get the relations, 
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These are the final expressions for the transverse and lon- 
gitudinal conductivities, respectively, which offer a very 
convenient method for numerical calculation of the con- 
ductivities using Kubo formalism. 

Throughout our study we choose the units where c = 
e = h = 1 and measure the energy scale in units of t. 



III. NUMERICAL RESULTS AND DISCUSSION 



A. Energy Spectrum 



both for the ordered as well as the disordered cases. The 
vertical lines correspond to the locations of the energy 
eigenvalues. For the ordered case we use black lines, 
while the blue lines are used for the disordered case, and, 
these two colored lines are superimposed in each plot to 
reveal the effect of disorder clearly. In the absence of 
disorder we get very sharp lines associated with the en- 
ergy eigenvalues of the lattice and all these energy levels 
are highly degenerate. The degeneracy factor of the en- 
ergy levels strongly depends on the specific choice of the 
magnetic field, i.e., on the value of Q. When Q becomes 
commensurate with N y , we get Q number of such sharp 
energy levels and each level becomes N/Q-iold degener- 
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FIG. 2: (Color online). Energy levels for a square lattice 
(30 x 30) when <f> is set at 0.1. The black lines represent the 
locations of the energy levels for the ordered (W = 0) case. 
On the other hand, the blue lines denote the positions of the 
energy levels for the disordered (W 7^ 0) case, where (a), (b) 
and (c) correspond to W = 1, 2 and 3, respectively. 



To make this present communication a self-contained 
one let us begin with the energy band structure of our 
model quantum system. 

In Fig. [5] we present the energy levels for a finite size 
square lattice in presence of a transverse magnetic field 



ate. These levels are quite similar to the Landau levels 
of the continuum model of 2D electron gas. It is to be 
noted that unlike the equispaced Landau levels of the 
continuum model, wc do not have equispaced energy lev- 
els for finite sized squarelatticc. The levels are widely 
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separated near the two edges of the energy spectrum as 
observed from Fig. [5J The situation is somewhat inter- 
esting when impurities are introduced into the system. 
In presence of disorder the degeneracy of each Landau 
band gets lifted due to the impurity potentials and the 
sharp energy levels form quasi-bands keeping the total 
number of energy levels in each quasi-band identical to 
each other and this number is equal to the degeneracy of 
the band. As far example, if we choose N x = N y = 30 
and set the magnetic flux (f> at 0.1 (i.e., Q = 10), then 
10 Landau bands appear in the spectrum (Fig. [2]) and 
each band contains 90 energy levels. The width of each 
Landau band increases with the increase of strength of 
disorder W as can be observed from the spectra, and for 
large enough disorder strength the neighboring Landau 
bands start to overlap with each other. 



of the Landau bands, the Hall conductivity enhances ex- 
actly by the same amount e 2 //i, and accordingly, for v 
filled Landau bands, o xy becomes equal to ve 2 /h, which 
is precisely the characteristic feature of integer quan- 
tum Hall effect. It is important to note that in tradi- 
tional ballistic waveguides we also get quantized conduc- 
tance, but such a strange precise quantization, about one 
part in million, can never be acheived in ordinary ballis- 
tic conductors as the backscattering processes are not 
completely suppressed. The almost complete elimination 
of the backscattering processes can only be obtained in 
the quantum Hall regime which actually ensures the ex- 
tremely high precision in the quantized value of Hall con- 
ductivity. Recently, it has been reported that when the 
Fermi energy lies within a plateau region, there is practi- 
cally zero overlap between the current carrying states in 
the sample leading to the almost complete suppression of 



B. Transverse and Longitudinal Conductivities 

Let us now illustrate the characteristics features of lon- 
gitudinal and transverse conductivities and the related 
issues for a finite sized square lattice. In Fig. [3] we show 
the variation of transverse conductivity a xy (red line) as 
a function of Fermi energy Ep for a (60 x 60) square 
lattice considering <f) = 0.1, W = 1 and fc^T = 0.01 
The Landau levels (light green) are also superimposed 




FIG. 3: (Color online). Transverse conductivity (red curve) 
as a function of Fermi energy for a square lattice (60 x 60) 
considering (f> = 0.1, W = 1 and ksT — 0.01. The Landau 
bands (light green) are superimposed on it. 

on the spectrum. From the spectrum it is clear that 
the transverse/Hall conductivity <j xy increases in inte- 
ger steps with the rise of Fermi energy Ef and at these 
plateaux the Hall conductivity gets the value, ve 2 /h with 
a great accuracy, where the integer number v corresponds 
to the total number of filled Landau bands below the 
Fermi level. For a particular system size the total num- 
ber of available Hall plateaux strongly depends on the 
choice of magnetic flux (j>, since the specific choice of 4> 
fixes the total number of Landau bands as mentioned 
earlier. Very interestingly we observe that even for a 
square lattice whenever the Fermi energy crosses any one 
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FIG. 4: (Color online). Longitudinal conductivity (red curve) 
as a function of Fermi energy for a square lattice (60 x 60) 
for the same parameter values given in Fig. [3] The Landau 
bands (light green) are superimposed on it. 

the backscattering processes and thereby the momentum 
relaxation process^. 

As a result of almost entire suppression of the backscat- 
tering processes for Fermi energy lying in the plateau 
regions, we would get practically zero longitudinal re- 
sistance and hence vanishing longitudinal conductivity 
which directly follows from the conductivity tensor. The 
result is shown in Fig. U where we set the same parameter 
values as given in Fig.[3J The resistance of a conductor is 
associated with the rate at which electrons can relax their 
momentum. To loose the momentum an electron has to 
be scattered through the allowed energy eigenstates in 
the bulk of the conductor. This almost vanishes when- 
ever the Fermi energy lies between two Landau bands. 
The longitudinal resistance is non-vanishing only when 
the Fermi energy is within a Landau band where the 
backscattering process are present. 

Figure [3] clearly shows that the Hall plateaux get ex- 
tended well inside the Landau bands and the Hall con- 
ductivity rises only within certain central region of each 
Landau bands. The underlying physics behind it as fol- 
lows. In presence of disorder Landau bands get broad- 
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FIG. 5: (Color online). Probability amplitude (PA) at different lattice sites (N) of a square lattice (60 x 60) for the same 
parameter values as mentioned in Fig. [3] The top and bottom panels correspond to the results for the eigenstates selected 
from the left and right edges of the Landau band, respectively, while the middle one is for the eigenstates lie in the centre 
of the Landau band. The three different columns are associated with the three different Landau bands shown in Fig. [jj] To 
understand the nature of energy eigenstates more clearly, in each figure we present the results for two states those are chosen 
from the respective regions. 



ened and the natures of all eigenstates in a Landau band 
are not identical. The states those lie along the two edges 
of the Landau band are localized, while the states in the 
middle are almost extended. To reveal the nature of the 
energy eigenstates whether they are quasi-localized or ex- 
tended, in Fig. [5]we show the variation of probability am- 
plitude at different lattice sites of a square lattice with the 
same set of parameters as those taken in Fig. [3] The top 
and bottom panels correspond to the energy levels chosen 
from the left and right edges of the three Landau bands, 
respectively, while the middle panel refers to some of the 
energy levels near the centre of the Landau bands. Quite 
clearly we see that near the band edges the probability 
amplitude vanishes almost at all sites excepting a few, 
which implies that electrons are localized in these states. 
On the other hand, for the states well inside the Lan- 
dau bands we have finite probability amplitude on each 
site revealing the extended nature of the states. This be- 
havior is true for every Landau band. Therefore, for the 
energy levels which lie near the two edges of the Landau 
bands we have almost zero contribution to the transverse 
conductivity and it leads to the plateau regions in the 
Hall conductivity, and, the non-zero contribution comes 
only from the extended states near the center of the Lan- 



dau bands. We also observe that the length of the Hall 
plateaux get widened with the rise of impurity strength, 
but it does not affect the quantized nature of Hall con- 
ductance. The quantization of Hall conductance is very 
robust and really very surprising. This is due to the fact 
that when wc increase the strength of disorder in the 
sample, more states become localized from the edges of 
the Landau bands and the regions of extended states in 
the band centre becomes narrow. However, the contri- 
bution to Hall conductivity that gets lost due to these 
additional localized states is exactly compensated by the 
enhanced contribution from the extended states, and ac- 
cordingly, the quantized nature of Hall conductivity is 
remains unaltered. This phenomenon is known as the 
current compensation and it has already been discussed 
in the literature by using the continuum models. 

C. Effect of Magnetic Field and Temperature 

So far we have discussed the main aspects of integer 
quantum Hall effect in a square lattice, but for the sake 
of completeness now we focus our attention on the ef- 
fects of magnetic field and temperature on IQHE as this 
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phenomenon is highly sensitive to the interplay of these 
parameters. With the increase of disorder strength the 
Hall plateaux can be found to be destroyed one by one 
m a generic sequence^. 

In Fig. [6] we present the variation of Hall conductiv- 
ity of a square lattice (60 x 60) as a function of Fermi 
energy for three different values of magnetic flux </> reg- 
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FIG. 6: (Color online). Hall conductivity as a function of 
Fermi energy for a square lattice (60 x 60) with W = 1 and 
KbT = 0.1. The dark- violet, pink and blue lines correspond 
to Q = 5, 10 and 15, respectively. 




FIG. 7: (Color online). Hall conductivity as a function of 
Fermi energy for a square lattice (60 x 60) considering W = 1 
and (j> = 0.1. The blue, pink and dark-violet colors correspond 
to fcsT = 0.01, 0.05, and 0.1, respectively. 

ulatcd by the parameter Q since 4> = 1/Q- Here we set 
ksT = 0.1 and W = 1. The dark- violet, pink and blue 
curves correspond to Q = 5, 10 and 15, respectively, and 
these values of Q are such that N x is exactly divisible by 



them so that the translational invariance remains along 
the x-dircction. For a given value of Q, we have Q Lan- 
dau bands in the energy spectrum, and accordingly, for 
a certain energy range the number of Hall plateaux are 
different for various Q as can be seen clearly from Fig. [6] 
However, in all these cases the quantized nature of the 
Hall conductance is exactly maintained, i.e., a xy is inte- 
ger multiple of the factor e 2 /h with a great accuracy. 

The effect of temperature in IQHE is illustrated in 
Fig. where we show the results for three temperatures 
considering W = 1 and flux <j) ~ 0.1. For a finite disorder 
strength and weak magnetic field, the widths of the Hall 
plateaux decrease gradually with the rise of temperature 
and the conductance spectrum becomes quite messy. For 
a sufficiently large temperature the phenomenon of con- 
ductance quantization almost disappears and this phe- 
nomenon has also been verified experimentally by some 
group o 20 i 21 . When the Fermi level is placed anywhere 
within the localized regions, the thermally excited elec- 
trons can jump toward the extended region and con- 
tribute to the conductivity. Thus, there is a competition 
between these two issues as a resultthe quatization gets 
lost gradually. 

IV. CONCLUSION 

In conclusion, we have re-visited the phenomenon of 
integer quantum Hall effect in a two-dimensional square 
lattice within a non-interacting electron picture. The in- 
terplay between the magnetic filed and randomness has 
been discussed in detail. We have used a single-band 
nearest-neighbor tight-binding Hamiltonian to describe 
the model quantum system and numerically calculated 
the longitudinal and transverse conductivities by using 
Kubo formalism. The physical picture about the integer 
quantum Hall effect that has emerged from our present 
investigation based on the discrete lattice model is al- 
most the same as that obtained in the continuum model. 
Our approach could be much more suitable for further 
investigation of the Quantum Hall effect in topologically 
insulating materials, like, graphene flakes, kagome lat- 
tices, etc. This is our first step towards this direction. 
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